End group effect on electrical transport through individual molecules: A microscopic 

study 
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The effect on molecular transport due to chemical modification of the metal-molecule interface 
is investigated, using as an example the prototypical molecular device formed by attaching a p- 
disubstituted benzene molecule onto two gold electrodes through chemically different end groups. 
Using a first-principles based self-consistent matrix Green's function method, we find that depending 
on the end group, transport through the molecule can be mediated by either near-resonant-tunneling 
or off-resonant-tunneling and the conductance of the molecule varies over more than two orders of 
magnitude. Despite the symmetric device structure of all the molecules studied, the applied bias 
voltage can be dropped either equally between the two metal-molecule contacts or mostly across 
the source (electron-injecting) contact depending on the potential landscape across the molecular 
junction at equilibrium. 

PACS numbers: 85.65.+h,73.63.-b,73.40.-c 



The continuing development of molecular electronics 
represents the convergence of the trend of device minia- 
turization and the growing expertise on single-molecule 
manipulation through scanning nanoprobcs, mechanical 
break junctions and supramolecular self-assembly tech- 
niques. 1 In contrast to carbon nanotubes and other in- 
organic semiconductor nanowires with uniform lattice 
structure^ devices based on individual organic molecules 
often require attaching appropriate end groups chemi- 
cally different from the molecule core in order to estab- 
lish stable contact to the metallic electrodesi The intro- 
duction of end groups into the molecular structure has 
two immediate consequences: (1) it introduces molecular 
states that are end group-based; (2) it modifies the metal- 
molecule interaction through the metal-end group bond. 
Single-molecule devices can therefore be considered as 
atomic-scale heterostructured devices, ^»» where the het- 
erostructure can be introduced at the metal-molecule in- 
terface through the end groups or in the molecule core 
through appropriate molecular design techniques. The 
purpose of this paper is thus to elucidate the prospect 
of interface "engineering" of molecular transport due to 
chemical modification of the metal-molecule interface, 
through detailed microscopic study of selected single- 
molecule devices. 

The interface "engineering" considered here is achieved 
through modifying the valence structure of the end 
group. The devices we consider are formed by attach- 
ing the p-disubstituted benzene molecule onto two semi- 
infinite gold electrodes through oxygen (O) and fluorine 
(F) end atoms and the isocyanide (C-N) end group. The 
device structure is shown schematically in Fig. In the 
case of C-N end group, the end atom can be either car- 
bon (C) or nitrogen (N). Since the hydroxyl substitucnt 
(O-H) may not deprotonate on contact with the metallic 
electrodes, we consider both O end atom and O-H end 
group in this work. For clarity in notation, we will de- 
note the five molecules following their structure as 0$0, 
HO$OH, F<f>F, CN$NC and NC$CN ($ stands for ben- 



zene ring) respectively. The molecules chosen are thus 
among the smallest possible where metal-molecule inter- 
action should have a strong effect on electrical transport, 
but are still representative of current experimental ef- 
forts i 

As the end group changes, both the energy and the 
charge distribution (symmetry) of the frontier molecu- 
lar orbitals (occupied and unoccupied molecular orbitals 
energetically closest to the metal Fermi- level, —5.31 
eV for gold(lll) electrode) change, depending on the 
chemical difference between the end group and the core 
atoms. For the NC$CN and CN$NC molecules, both the 
highest-occupied-molecular-orbital (HOMO) and lowest- 
unoccupied-molecular-orbital (LUMO) are delocalized 
through the entire molecule (central benzene as well as 
end group) and involve mainly the orbitals of carbon 
and nitrogen. Switching to O and F atoms, the HOMOs 
remain delocalized through the entire molecule, but the 
LUMOs become localized within the molecular core (cen- 
tral benzene). The HOMO-LUMO gap increases as we 
move from O to F (this is expected since these electron- 
withdrawing substituents stabilize the occupied levels). 
The frontier molecular states of the 0<f>0 molecule (spin- 
singlet) show similar features to the phenyl dithiolate 
(PDT or S<&S in the present notation) molecula&i due 
to the identical valence structure of the O and S atoms. 
Adding the H end atom increases the HOMO-LUMO gap 
of the HO$OH molecule. Although the HOMO of the 
HO$OH molecule remains delocalized over the oxygen 
and carbon atoms, it has little weight on the H ends 
which effectively reduces its coupling with the metal sur- 
face states. 

Due to differences in the energy and charge distribution 
of the frontier molecular states, the above molecules show 
distinctly different behavior in the molecular response 
to the perturbation induced by the metal-molecule in- 
teraction and the applied bias voltage. This is in- 
vestigated using a first-principles based self-consistent 
matrix Green's function theory which combines the 
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Non-Equilibrium Green's Function theory of quantum 
transport with an effective single-particle description of 
molecular junction electronic structure using density- 
functional theory (DFT)M& The theoretical models used 
have been discussed extensively in our previous workj&4 
where we presented the general theoretical approach 
for modeling transport through single-molecule devices 
within the coherent transport regime and identified two 
key device processes for understanding the transport 
characteristics of a two-terminal molecular device: the 
equilibrium energy-level lineup and the noncquilibrium 
charge/potential response to the applied bias. The same 
theoretical and modeling approaches are adopted here. 
To summarize, the gold electrodes are modeled as semi- 
infinite (111) single-crystals. Six nearest-neighbor gold 
atoms on each metal surface (twelve gold atoms over- 
all) are included into the "extended molecule" where the 
self-consistent calculation is performed. The rest of the 
electrodes (with the six atoms on each surface removed) 
are considered as infinite electron reservoirs, whose ef- 
fects are included as self-energy operators. The calcu- 
lation is performed using a modified version of Gaus- 
sian98 5 using the Becke-Perdew-Wang parameterization 
of density-functional theory— and appropriate pseudopo- 
tentials with corresponding optimized Gaussian basis 
sets^i For comparison with our previous work and other 
experimental and theoretical works^^ we keep the same 
adsorption geometry and metal surface-end group dis- 
tance for all the molecules studiedA^ The calculation is 
performed at room temperature. 

Following Ref. 0, we analyze the device physics of the 
metal-molecule-metal junction both at equilibrium and 
out of equilibrium. The most important quantity at 
equilibrium is the charge transfer between the molecule 
and the electrodes upon formation of the metal-molecule- 
metal junction, which may be understood qualitatively 
from the bonding configuration change across the metal- 
molecule interface and analyzed through the symme- 
try and energy of the frontier molecular orbitals. The 
transferred charge and the induced change in the elec- 
trostatic potential are obtained by taking the difference 
between the self-consistent charge/potential distribution 
in the equilibrium metal-molecule-metal junction and 
the charge/potential distribution in the isolated molecule 
plus the bare bimetallic junction. For the CN<I>NC and 
NC^CN molecules, electron densities increase on the the 
end atom P^ orbitals (Fig. 0). For the 04>0 molecule, 
the overall spatial distribution of the transferred charge 
is similar to that of the PDT molecule (with electron 
density increasing on the end atom P^ orbitals but de- 
creasing in the P x orbitals)^ but the magnitude of the 
charge transfer is larger due to the larger electronegativ- 
ity of the oxygen atom which leads to a larger barrier 
for electron injection from the metal into the molecule 
(Fig. 0) . Keeping the H atom reduces significantly the 
amount of charge transferred to the HO<i>OH molecule 
due to the saturated H-0 bond (Fig. Since most of 
the charge transfer is to the end H and O atoms, this 



leads to a large potential barrier at the H-0 bonding re- 
gion (Fig.[3J). The F$F molecule shows different behav- 
ior compared to the other molecules. The lowest unoc- 
cupied state is core-based, and is located about 3.8(eV) 
above the metal Fermi-level. So the bonding configura- 
tion change is mainly realized by the charge redistribu- 
tion among the occupied molecular states, leading to an 
electron density decrease in the fluorine-carbon a bond 
and an electron density increase in the fluorine-surface 
bond. The magnitude of the charge transfer is much 
smaller than the 0<&0 molecule and gives rise to an elec- 
trostatic potential well instead of electrostatic potential 
barrier for electron injection into the molecule (Fig.[3J). 

The charge transfer processes and the resulting change 
in the electrostatic potential determine the energy-level 
lineup relative to the metal Fermi-level for the equilib- 
rium molecular junction as well as low-bias conductance, 
which is obtained from the electron transmission charac- 
teristics at zero bias and calculated using the method 
of Refs. 13 and H The results are illustrated in Fig. 
^ For comparison we have also shown the transmis- 
sion characteristics of the PDT molecule. For both the 
CN$NC and NC$CN molecules, the metal Fermi-level 
Ep is in near-resonance with the LUMO state, which is 
delocalized over the entire molecule. This leads to broad 
transmission spectrum around Ep and large zero-bias 
conductance (close to the conductance quantum) due to 
near-resonant-tunneling through LUMO. The difference 
in the transmission spectrum above Ep between the two 
molecules is due to the difference in the energy-level spac- 
ing and broadening of the unoccupied states (the unoccu- 
pied states of CN$CN are more closely spaced and more 
broadened than those of NC$CN). The HOMO-LUMO 
gap of NC$CN is also larger than that of CN$CN, push- 
ing the HOMO-mediated transmission peak further be- 
low Ep. For the 0$0 and F$F molecules where the 
end group atoms are more electron-rich, transmission at 
the metal Fermi-level is due to tunneling through the 
HOMO-LUMO gap. Although transmission through the 
middle of the gap is small, the HOMO of the 0$0 
molecule lines up much closer to the metal Fermi-level 
than PDT, leading to a much larger conductance. Keep- 
ing the end H atoms (HO<f>OH) pushes both peaks cor- 
responding to HOMO and LUMO down, with Ep ly- 
ing slightly closer to HOMO than LUMO, which reduces 
the conductance of the HO<I>OH junction by a factor of 
1=3 100. Finally for the F$F molecule where the end atom 
has the largest number of valence electrons, the LUMO 
instead lines up closer to the metal Fermi-level than the 
HOMO. Since the LUMO state of F$F is localized on the 
molecule core and the HOMO state is pushed far below 
the metal Fermi-level (beyond the energy range plotted 
in Fig. 2J), the zero-bias conductance is very small. As 
the end group changes, the conductance of the molecu- 
lar junction varies by more than two orders of magnitude, 
from the maximum of 45.8 (nS) for the CN<I>NC molecule 
to 0.18 (/zS) for F$F molecule. 

For the device out of equilibrium, the central quan- 
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tity is the bias-induced charge redistribution within the 
molecular junction and the spatial distribution of the 
voltage drop, which determine the current-voltage (I- 
V), differential conductance- voltage (G-V) characteris- 
tics and the bias-induced modification of molecular states 
(the static Stark effect). The I-V and G-V characteris- 
tics are calculated self-consistently at each bias voltage 
and shown in Fig. For the CN$NC and NC$CN 
molecules, the conductance peaks near zero-bias are 
due to near-resonant-tunneling through the LUMO. For 
the 0<I>0 molecule, it is due to near-resonant-tunncling 
through the HOMO. Within the bias range considered 
here, both the F<i>F molecule and the HO<I>OH molecule 
arc insulating due to the small tunneling probability 
through the HOMO-LUMO gap, leading to nearly-linear 
(Ohmic) I-V characteristics from — 1(V) to +1(V). By 
examining the transmission characteristics as a function 
of the applied bias (not shown here) , we find that the in- 
crease of conductance towards ±2(V) is due to the closer 
alignment of the metal Fermi-levels with the LUMO 
for the F$F molecule and the HOMO for the HO$OH 
molecule. 

Although the device structures considered here are all 
symmetric with respect to the exchange of the source 
and drain electrodes, the voltage drop across the metal- 
molecule-metal junction at finite bias voltage shows dif- 
ferent behavior depending on the potential landscape 
throughout the equilibrium molecular junction (the volt- 
age drop is obtained by evaluating the difference be- 
tween the electrostatic potential at finite and zero bi- 
ases, which obeys the boundary condition of approaching 
±V/2 inside the electrodes). For the CN$NC, NC$CN 
and HO$OH molecules where there is a narrow po- 
tential barrier at the end group region but not in the 
molecule core at equilibrium, the voltage drops (nearly) 
equally between the two metal-molecule interfaces. But 
for the 0<I>0 molecule where the potential barrier for 
electron injection is much larger and extends through- 
out the molecule (end group plus molecule core) region 
at equilibrium, most of the voltage drop occurs at the 
interface between the molecule and the source (electron- 
injecting) electrode, as illustrated in Fig. EJfor bias volt- 
age of 2(V). This is readily analyzed using the concept 
of resistivity dipole for nonequilibrium transport systems 
in the presence of current flow^*^ where charge dipoles 
and correspondingly local electric fields develop in the 
vicinity of scattering centers (potential barriers here) to 
ensure the current continuity throughout the conduction 
system. For the 0$0 molecule with a large and thick 



potential barrier for electron injection, a large electrical 
field is required across the molecule-right electrode in- 
terface at positive bias voltage, which is sustained by a 
large electron deficit on the molecule side and leads to a 
large amount of the voltage drop across the correspond- 
ing interface (similar situation occurs at the molecule- 
left electrode interface at negative bias). In the case of 
the CN$NC, NC$CN and HO$OH molecules, only a 
relatively small electric field is needed for electrons to 
penetrate a narrow barrier at the molecule-source elec- 
trode interface, so the electron redistributions around 
the two metal-molecule interfaces have similar magnitude 
and voltage drops nearly equally across the two inter- 
faces. For the F$F molecule where there is an electro- 
static potential well in the equilibrium junction, the volt- 
age also drops equally between the two metal-molecule 
interfaces due to the absence of barrier. In addition, 
the difference in the potential response of the molecules 
to the applied bias voltage leads to different behavior 
in the bias-induced shift of molecular levels. For the 
CN$NC, NC$CN, HO$OH and F$F molecules, the en- 
ergy levels of HOMO and LUMO are nearly constant 
within the bias voltage range considered here. But for 
the 0<i>0 molecule where the voltage drops mostly across 
the molecule-source electrode interface, the shift of the 
HOMO and LUMO levels follows the Fermi-level of the 
drain (electron-extracting) electrodes. 

Discussion and Conclusion Most current experimen- 
tal works on molecular electronics have focused on 
molecules attached to the gold electrodes through sul- 
fur end groupfii due to the convenient thiol-gold self- 
assembly scheme. In this work, we have found theo- 
retically that depending on the end groups used, the 
molecules with identical molecule cores can show dis- 
tinctly different behavior in their response to the per- 
turbation induced by the metal-molecule interaction and 
the applied bias voltage, leading to widely different trans- 
port characteristics. Since one of the major advantages 
of molecular electronics is the potential to build devices 
with the desired properties from the bottom up, it is 
important to explore the feasibility of such device engi- 
neering through molecular design of interfaces. Although 
we are not aware of any existing single-molecule measure- 
ments that test the end group effects as reported here^ 
such tests should be plausible, in particular using the 
mechanical break junction techniques^ 
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FIG. 1: (Color online) Schematic illustration of the metal- 
molecule-metal junction. Six gold atoms on each metal sur- 
face are included into the "extended molecule" where the self- 
consistent calculation is performed. The molecule can be at- 
tached onto the electrodes through chemically-different end 
groups. 
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FIG. 2: (Color online) Charge transfer (in unit of l/(a.u.) 2 ) 
and cross-sectional view of the electrostatic potential change 
upon the formation of the gold-CN$NC-gold (left figure) and 
gold-0$0-gold (right figure) junctions. Here the spatial dis- 
tribution of the transferred electrons is plotted as a function 
of position in the xy-plane (denned by the benzene ring) after 
integrating over the z-axis. Also shown is the projection of 
the molecule onto the xy-plane. 




FIG. 3: (Color online) Charge transfer (in unit of l/(a.u.) 2 ) 
and cross-sectional view of the electrostatic potential change 
upon the formation of the gold-F<E>F-gold (left figure) and 
gold-HO<£>OH-gold (right figure) junctions. Here the spatial 
distribution of the transferred electrons is plotted as a func- 
tion of position in the xy-plane (defined by the benzene ring) 
after integrating over the z-axis. Also shown is the projection 
of the molecule onto the xy-plane. 
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FIG. 4: Transmission versus energy characteristics of the 
equilibrium molecular junction (at zero bias). The horizontal 
lines show the location of the metal Fermi-level. 




FIG. 5: Current-voltage (upper figure) and differential 
conductance-voltage (lower figure) characteristics of the five 
molecules studied in this work. 





FIG. 6: (Color online) Charge redistribution and cross- 
sectional view of the voltage drop across the gold-CN$NC- 
gold (left figure) and the gold-0«J>0-gold (right figure) junc- 
tion at bias voltage of 2(V). The transferred charge is ob- 
tained by integrating the difference in the electron density at 
finite and zero biases and plotted as a function of position in 
the xy-plane. 



